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Motivated by the recent discovery of superconductivity in the iron-based ladder compound 
BaFe 2 S 3 under high pressure, we derive low-energy effective Hamiltonians from first principles. 

We show that the complex band structure around the Fermi level is represented only by the Fe 3d xz 
(mixed with 3 d xy ) and 3d x 2_ y 2 orbitals. The characteristic band degeneracy allows us to construct 
a four-band model with the band unfolding approach. We also estimate the interaction parameters 
and show that the system is more correlated than the 1111 family of iron-based superconductors. 
Provided the superconductivity is mediated by spin fluctuations, the 3d;cz-like band plays an essen¬ 
tial role, and the gap function changes its sign between the Fermi surface around the T point and 
that around the Brillouin-zone boundary. 


I. INTRODUCTION 

Since the discovery of superconductivity in fluorine 
doped LaFeAsOi, a variety of iron pnictides and chalco- 
genides have been found to exhibit superconductivity 
with high transition temperatures (T c ). In these com¬ 
pounds, Fe ions commonly form a two-dimensional (2D) 
network. This fact raises an intriguing and fundamental 
question whether the square network is essential for the 
high T c superconductivity and what happens in different 
geometries. In the case of the cuprates, it was discovered 
that the ladder compound (Sr,Ca)i 4 Cu 2404 i becomes a 
superconductor under pressure ~ 3 GPa.— This exper¬ 
imental observation has stimulated various theoretical 
studies on superconductivity in quasi-ID systems. In¬ 
terestingly, it has been recently found that a ladder com¬ 
pound BaFe 2 S 3 becomes a superconductor under high 
pressure ~ lOGPa.— In the phase diagram, the super¬ 
conducting phase resides next to a magnetic insulating 
phase, in which the spin correlation is antiferromagnetic 
along the ladder and ferromagnetic along the rung. The 
maximum T c is as high as 14 K. The superconducting 
transition in BaFe 2 S 3 is of great interest, since we may 
have a chance to pin down the origin/mechanism of the 
high T c superconductivity in iron-based superconductors 
(FeSC) by investigating the commonalities and differ¬ 
ences between the 2D and quasi-ID systems. 

Recently, motivated by the experimental works on su¬ 
perconductivity in the single-layer potassium-doped iron 
selenide (110) film (which can be viewed as a weakly 
coupled ladder system)^ or ladder compounds such as 
BaFe 2 Se 3 ^r— and CsFe 2 Se;?ii (which do not exhibit su¬ 
perconductivity), a variety of theoretical studies have 
been reportedHowever, studies for BaFe 2 S 3 based 
on ab initio calculation are yet to be performed. In 
this study, we derive low-energy effective Hamiltonians 
for BaFe 2 S 3 from first principles, focusing on the Fe 3d 
bands around the Fermi level (Ef)- We find that a d xz - 
like orbital, a linear combination of d xz and d xy , forms 


two Fermi (electron) pockets around k z =0 and a d x 2 _ y 2 - 
like orbital forms a Fermi (hole) pocket around k z =ir. 
The effective energy bands obtained with the two or¬ 
bital models are further depicted in the unfolded Bril- 
louin zone (BZ), in which one of the electron pockets is 
placed around the BZ boundary of k' z =n. Since the mag¬ 
netic instability is strong at q' z ~ n in the extended BZ, 
the d X2 -like orbital should be active for spin-fluctuation- 
mediated superconductivity. We also estimate the values 
of interaction parameters in the effective Hamiltonian, 
such as the Hubbard U and Hund’s coupling J. We show 
that BaFe 2 S 3 is more strongly correlated than the 1111 
compounds. 


II. CRYSTAL STRUCTURE 


In Fig. [TJ we show the crystal structure of BaFe 2 S 3 . 
We see that Fe atoms form two-leg ladders running along 
the c-axis. In the following calculation for ambient pres¬ 
sure, we used the lattice constants a, b , c and the atomic 
positions of Ba, Fe and S reported in Ref. [l9|. Namely, 
the lattice constant a, b, and c are 8.78, 11.23 and 5.29 
A, respectively. The space group is Cmcm, and the 
atomic positions of Ba(4c), Fe(8e), S(4c) and S(8g) are 
(0.0, 0.686, 0.25), (0.154, 0.0, 0.0), (0.0, 0.116, 0.25) and 
(0.208, 0.378, 0.25), respectively. In the phase diagram of 
temperature (T) and pressure (P), the superconducting 
phase has a dome-like shape and T c is highest around 
P=12.4 GPa. For P=12.4 GPa, a, b and c shrink to 
96.0%, 92.0% and 96.6% of those at ambient pressure, 
respectivel y 19 ’ 20 . Since the atomic configuration under 
pressure is yet to be reported, in the present calculation, 
we just change the lattice constants. 
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FIG. 1: (color online): Crystal structure of BaFe 2 S 3 . 


III. BAND STRUCTURE, FERMI SURFACE 
AND PROJECTED DENSITY OF STATES 


(a) P = 0 GPa 



(b) P = 12.36 GPa 
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We performed density functional calculation for the 
experimental crystal structure^ with the QUANTUM 
ESPRESSO package * 22 Here we employ the exchange corre¬ 
lation functional by Perdew et al 2^. The wave functions 
are expanded by plane waves up to a cutoff energy of 40 
Ry, and an 8 x 8 x 8 fc-mesh in the first Brillouin zone is 
used. 

Based on the band calculation, we further con¬ 
structed Wannier functions for these bands using the 
WANNIER90 package^. The resulting spreads of the 
d^ z 2 _ r 2 : d xz ,dy Z ,d x 2 _ y 2 , and d xy orbitals are 2.35, 3.51, 
3.13, 2.40 and 3.14 A 2 for ambient pressure, and 2.91, 
4.01, 3.35, 2.73 and 3.43 A 2 for P = 12.4 GPa. If we 
compare the Wannier spreads of the 2D FeSC (Table III 
in Ref. l25l ). we see that at P =0 the Wannier functions 
are as localized as those of LiFeAs, but at P=12.4 GPa 
they are more delocalized as those of BaFe 2 As 2 - 

In Fig. [2j we show the band dispersion obtained by 
the Wannier interpolation (blue dashed curves) with that 
from the density functional calculation (red solid curves). 
Since the unit cell contains four Fe atoms, there are 
twenty Fe 3d bands around the Fermi level. We see that 
the band width of the Fe 3d states for P=12.4 GPa is 
larger than that for P=0 by ~ 25%. We list the resulting 
hopping integrals in the Supplemental Material^ As is 
expected from the crystal structure shown in Fig. [T| the 
band structure is dispersive along the k z axis (from the 
r to Z point). The transfer integral is largest between 
the nearest neighbor d 3z 2_ r 2 orbitals along the ladder 
(the c-axis), and its amplitude is ~ 0.56 eV. For P=12.4 
GPa, it becomes 0.64 eV. These values are larger than 
the nearest-neighbor transfer integrals for the 2D FeSC 
(see Tables IV-VII in Ref. [25l ). On the other hand, the 
inter-ladder transfer integrals are small, so that the band 
structure of BaFe 2 S 3 is quasi-one dimensional along the 
k z axis. Nevertheless, since the dominant transfer in¬ 
tegrals along the c axis are large, the band width of 
BaFe 2 S 3 turns out to be a similar value to that of 2D 


FIG. 2: (color online): Band structure for ambient pressure 
(a) and P=12.4 GPa (b). The original (Wannier interpolated) 
band dispersion is shown by red solid (blue dashed) curves. 
In the inset, we show the Brillouin zone. 



FIG. 3: (color online): Fermi pockets a and {3 (a) and 7 (b) 
of BaFe 2 S 3 for P=0. 

FeSC (see Fig. 4 in Ref. I25I ). As for the onsite energies 
of d 3z 2 _ r 2 ,d xz ,d yz ,d x 2 _ y 2 , and d xy , they are 7.44, 8.00, 
8.03, 7.69, and 7.95 eV (where E F = 8.19 eV) for P=0, 
and 9.37, 10.02, 9.98, 9.62 and 9.89 eV (where E F = 10.22 
eV) for P=12.4 GPa. The size of the crystal field split¬ 
ting is similar to that of the 2D FeSC. 

In Figj3l we show the Fermi surface for P=0. There 
are two electron pockets around k z = 0 , hereafter called 
a and /3, and one hole pocket around k z = 7 r, called 
7 . For P=12.4 GPa, the /? pocket significantly shrinks, 
while the a pocket is quite robust against the external 
pressure. The energy band lying just below E F {=0 eV 
hereafter) around the T points at P=0 crosses the Fermi 
level under pressure, resulting in another tiny pocket for 
P=12.4 GPa. 

We plot the projected density of states for the 3d or- 
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FIG. 4: (color online): Projected density of states for P= 0. 
While the amplitude of d yz and d 3z 2_ r 2 around Ef is not 
large, the contributions from d x 2 _ y 2 , d xz and d xy are signifi¬ 
cant. 



FIG. 5: (color online): (a) Band dispersion of the effective 
two-orbital (eight-band) model for P=0 (blue solid curves). 
The original ab initio band structure is shown by dotted 
curves, (b) Band dispersion in the extended (unfolded) BZ. 


bitals in Fig. SI We see that d 3z 2_ r 2 does not contribute 
to the low-energy states. On the other hand, the Fermi 
pockets ct, /3, and 7 are mainly composed of d xy , d xz and 
d x 2 _ y 2 , respectively. The contribution of d yz to these 
pockets is subdominant. As for the electron filling n 
of d 3z 2 _ r 2 1 d xzi d yz ,d x 2 _ y 2 , and d xy , they are 0.61, 0.54, 
0.56, 0.76, and 0.50 for P=0, and 0.60, 0.51, 0.62, 0.76, 
and 0.53 for P=12.4 GPa. Since n of the d x 2_ y 2 orbital 
is quite high, d x 2_ y 2 should be inactive magnetically. It 
is interesting to note that n of d 3z 2_ r 2 in the 2D FeSC is 
~ 0.75 (see Table II in Ref. [ 2 ^), and the orbital extends 
in the direction perpendicular to the Fe plane as d x 2_ y 2 
in BaFe 2 S 3 does. So far, several studies have pointed out 
that d 3z 2_ r 2 in the 2D FeSC does not play a crucial role 
in the magnetism and superconductivity—r— 


In the 2D FeSC, both the hole and electron pockets 
consist of more than two orbitals. Here let us look at 
the orbital character of the Fermi surface of BaFe 2 S 3 . 
In Fig® we plot the band dispersion of the four-band 
model. The contribution of the each orbital is repre¬ 
sented by the width of the curves. We can see that w xz 
forms electron pockets around T and Y' (the Fermi pock¬ 
ets a and f3 in the original BZ). On the other hand, the 
hole pocket sitting between T and T' (the 7 pocket in the 
original BZ) consists of w xz and w x 2 _ y 2 . The far (near) 
side of the pocket from T has a dominant w x 2 _ y 2 - (w xz ~) 
character. 


V. LINDHARD FUNCTION 


IV. TWO-ORBITAL MODEL AND BZ 
UNFOLDING 

We now further try to simplify the five-orbital (twenty- 
band) model. We have found that the energy dispersion 
around Ep can be represented by a two-orbital (eight- 
band) model in which the Wannier functions are mainly 
composed of two orbitals, i.e. d x 2_ y 2 and d xz hybridized 
with d xy . We hereafter call the Wannier functions w x 2_ y 2 
and w xz . To construct the Wannier functions, the inner 
(outer) window is set to be [-0.3 eV,0.2 eV] ([-1.2eV,1.5 
eV]) with respect to Ep — In Fig. ®a), we compare the 
original ab initio bands and the Wannier interpolated 
bands for ambient pressure. We obtained a similar re¬ 
sult for P=12.4 GPa (not shown). We list the resulting 
hopping integrals for the two-orbital model in the Sup¬ 
plemental Material— 

Interestingly, we can further simplify the current 
model by unfolding the BZ, as was done for the 1111 
compound— Indeed, if we introduce a local gauge trans¬ 
formation for one of the two orbitals to change its sign, 
we can expand the band dispersion from r to Z (see 
Fig®b))— Then the model is simplified to a four-band 
model. 


In the 2D FeSC, the Lindhard function x°(q) has a 
peak around (7r,0,0) and ( 0 , 7 r, 0 ) in the extended BZ^i, 
which is compatible with the stripe-type antiferromag¬ 
netic order. This fact has caused hot debates on which 
of the localized spin picture or the itinerant picture is 
more appropriate for describing the magnetic properties 
of 2D FeSC — — Thus it is interesting to see whether the 
peak position of x°(q) for BaFe 2 S 3 agrees with the Bragg 
peak observed in the experiment. 

In Fig® we plot x?j(q) f° r the five-orbital model in 
the extended Brillouin zone. Here i and j (=1,2) specify 
the Fe atom in the unit cell. We see that the correlation 
along the ladder (i = j) is stronger than that along the 
rung i ^ j, and x°(q) does not depend on q x and q y 
significantly. As for the q z dependence, x°(q) has larger 
values in the plane of q z = it than q' z = 0 , due to the 
nesting between the a and (3 pockets in Fig® Note that 
the (3 pocket moves to the plane of q' z = tt in the ex¬ 
tended Brillouin zone. For P=12.4 GPa, although the /3 
pocket shrinks, we confirmed that the contribution from 
the particle-hole scattering between a and is still dom¬ 
inant. This q' z dependence of x°(q) is consistent with the 
experiment in that the intra-orbital spin correlation is an¬ 
tiferromagnetic along the ladder. As for the correlation 
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FIG. 6: (color online): Band dispersion of the four-band 
model for P= 0. The weight of the each orbitals is orbital 
are represented by the width of the curve. 


VI. SUPERCONDUCTIVITY 


Since the superconducting phase resides next to the 
antiferromagnetic phase, let us here assume that the su¬ 
perconductivity is mediated by spin fluctuations. Then 
the pairing interaction is strong at q' z = 7r, and favors 
a pairing gap function having different signs between 
the Fermi pockets a and /3. The active band for su¬ 
perconductivity is w xz in this scenario, because both a 
and /3 are made from w xz . On the other hand, w x 2_ y 2 
makes a pocket around q' z = ±7r/2, but should be pas¬ 
sive for the superconductivity. In the 2D FeSC, it has 
been proposed that d xy is important to understand the 
material dependence of the superconducting transition 
temperature 42r— It is interesting to note that both w xz 
in BaFe 2 S 3 and d xy in the 2D FeSC extend in the direc¬ 
tion of pnictogen/chalcogen sites, and almost half-filled 
in the five-orbital model. 

For 2D FeSC, the orbital fluctuation mechanism has 
been extensively studied^ - — While it is an interesting 
future problem to study whether orbital fluctuations me¬ 
diate superconductivity in BaFe 2 S 3 , there is a notable 
difference between the quasi-ID and 2D systems. In the 
case of the 2D FeSC, d yz and d xz are essential for the 
orbital-fluctuation mechanism. While these orbitals cor¬ 


respond to d xy and d yz in BaFe 2 S 3 , these are irrelevant 
in the two-orbital model. 


.cr 
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FIG. 7: (color online): Wave number dependence of the Lind- 
hard function for the five-orbital model at P= 0. 


along the rung, it is ferromagnetic, since ferromagneti- 
cally coupled spin configuration S± + S 2 (where S', is the 
spin at the i-th Fe atom in the unit cell) has stronger 
correlation than Si — S 2 (note that X 12 > 0 at q = Y’). 
This is also consistent with the experiment. 

On the other hand, for the inter-ladder correlation, 
X°(q) is peaked at Y', while the Bragg peak in experi¬ 
ment is observed to be located at 0.5bi+0.5b 2 (bi and 
b 2 are the reciprocal primitive vectors) * This result sug¬ 
gests that the inter-ladder correlation can not be ex¬ 
plained simply in terms of Fermi surface nesting. 


VII. ELECTRON CORRELATIONS 


Let us move on to the electron correlations in BaFe 2 S 3 . 
By means of the constrained random phase approxima¬ 
tion (cRPA)£i, we estimate the values of the Hubbard U 
and the Hund’s coupling J in the five-orbital model. We 
used the density response code of El k 52 : 53 . In the calcula¬ 
tion of the charge susceptibility, we took 100 unoccupied 
bands and 4x4x4 k and q meshes. The double Fourier 
transform of the charge susceptibility was done with the 
cutoff of |G + q|=2,3,4,5 and 6 (1/a.u.) with G being 
the reciprocal vector, from which we estimated the values 
extrapolated to the limit of |G + q| —► 00. In Table Q] we 
list the resulting values of U and J for P = 0. If we com¬ 
pare these values (and the ratio J/U, which is important 
to measure the electron correlations in 2D FeSC^i) with 
those of the 2D FeS C 25 ’ 55 , we see that they are as strong 
as those of LiFeAs. The pressure dependence of U and J 
is not so strong and they become smaller only by ~ 6-7% 
under the pressures up to P=12.4 Pa. These results sug¬ 
gest that BaFe 2 S 3 is more strongly correlated than the 
1111 compounds : 25 ’ 56 This observation is consistent with 
the fact that the ordered moment of BaFe 2 S 3 is about 1.3 
fiB, which is larger than that of LaFeAsO 56 ' 57 . 
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TABLE I: Hubbard U and Hund’s J in the five-orbital model 
for P= 0 GPa. 


U [eV] 

xy 

yz 

3 z z -r 2 

XZ 

2 2" 

x -y 

xy 

3.35 

2.42 

2.34 

2.33 

3.16 

yz 

2.42 

3.26 

2.89 

2.30 

2.54 

3z 2 -r 2 

2.34 

2.89 

3.79 

2.77 

2.44 

xz 

2.33 

2.30 

2.77 

2.99 

2.43 

x 2 -y 2 

3.16 

2.54 

2.44 

2.43 

3.86 

J [eV] 

xy 

yz 

3 z 2 -r 2 

XZ 

2 2 

x -y 

xy 

- 

0.50 

0.67 

0.48 

0.22 

yz 

0.50 

- 

0.32 

0.46 

0.55 

o 2 2 

3z -r 

0.67 

0.32 

- 

0.30 

0.71 

xz 

0.48 

0.46 

0.30 

- 

0.51 

2 2 
x -y 

0.22 

0.55 

0.71 

0.51 

- 


VIII. CONCLUSION 

We have studied the electronic structure of BaFe 2 S 3 
from first principles. Detailed analysis of the full Fe 
orbital model revealed that the Fermi surfaces are rep¬ 


resented only by the two Fe orbitals, i.e. w xz and 
w x 2 _ y 2 . Provided that the superconductivity is medi¬ 
ated by spin fluctuations, w xz , which corresponds to 
d xy in the 2D FeSC, is essential for superconductivity. 
In fact, they have the following common features: (i) 
they extend in the direction of pnictogen/chalcogen sites, 
(ii) they are nearly half-filled, (iii) they form discon¬ 
nected Fermi pockets around T and near the boundary 
of the BZ, between which pairing interaction mediated 
by stripe-type antiferromagnetic spin fluctuation works 
effectively^ Concerning the electron correlations, it is 
stronger than those of LaFeAsO or LaFePO. 
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